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SPSECT OP VARIABLE THERl'^lAL CONDUCTIVITY AND HEAT 
TRANSFER COEPPICIENT ON PBRPORMANCB OP PINS 

This work presents a theoretical investigation into 
the effect of variable thermal conauctivity and heat 
transfer coefficient on the excess temperature distribu- 
tions, temperature gradients and the efficiencies of fins 
of three different shapes, A general differential equation 
governing the heat transfer by conduction and convection 
in fins is derived. A parabolic expression for the vari- 
ation of the temperature dependent thermal conductivity and 
an experimentally verified power series type of an expression 
for the heat transfer coefficient have been assumed. 

Results have been compared with those of the simple case, in 
which both the thermal conductivity and the heat transfer 
coefficient are held constant. The results reveal that 
variable thermal properties exercise a significant effect 
on the -performance of the fins considered. 



CHAPTER 1 


IHTROEUCTION 


1.1 INTRODUCTION 

In a conventional heat exchanger beat is transferred 
from one fluid to another through a metallic wall and, 
other things being equal, the rate of heat flow is directly 
proportional to the extent of the wall surface and to the 
temperature difference between one fluid and the adjacent 
surface. If thin strips of metal are attached to the basic 
surface, extending into one of the fluids, the total surface 
is thereby increased which results in an enhanced rate of 
heat flow. These attached heat conducting strips constitute 
what are generally termed ’’extended surfaces” or "fins”. 

Dins are of various types, the most common among them 
being (a) Rectangular , (2 ) Annular, and (3) Spines, 

The exact analysis of the heat transfer process in 
fins is a complicated problem. Analytical solutions 
have, however, been obtained after invoking a set of 
simplifying assumptions, essentially put forward by Murray [1 ] 
and G-ardner [2], which posit thermal properties namely 
the thermal conductivity and the heat transfer coefficient 
as constant. However both of these assumptions are open 
to serious questions. Temperature dependence, of thermal 
conductivity is sometimes not only strong but also bizarre [ 3 ]. 



There is also evidence suggesting a rapid increase in 
the magnitude of heat transfer coefficient towards the 
end of a fin and its variation from point to point along the 
length of the fin [4-7 ]. Against this background an attempt 
at a more rigorous analysis aimed at ascertaining the impact 
of these variable properties on fin performance seems to be 
in order, 

1.2 LITERATURE REVIStf 

Although the literature pertaining to the study of 
fins dates back to the work of I@:|^ehouse [s] in 1789, most 
of the work since has been concentrated on convecting-fins 
with constant thermal properties. It is only lately that 
the effect of variable thermal properties on the performance 
of fins, whether convecting, radiating, or both, has been 
studied and shown to be quite significant particularly when 
large temperature difference exists. Cobble [ 9 ] in 1963 
found the temperature distribution of a fin transferring 
heat by combined convection and radiation. He obtained 
the solution in terms of Jacobian elliptical functions 
with the aid of the Gregory-Newton forward interpolation 
formula, Stockman and Kramer [10 ] , in the same year, 
presented the effects of a linearly varying conductivity 
and anissivity on the heat transfer of purely radiating 
fins, using the Runga-Kutt a method. Hung and Appl [11] 
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analyzed fins vith simultaneous surface convection and 
radiation and presented the temperature distribution for 
variable thermal properties and heat generation. Kosarev 
and Nevrovskii [12] considered a purely convecting-fin 
with temperature dependent theimal conductivity aM developed 
the method of presenting the tip temperature ani the heat 
transfer rate. A regular perturbation solution for a 
straight convecting fin of constant thickness with temper- 
ature dependent conductivity has been presented by Aziz 
and Huq [13I More recently Razelos and Imre [14] found 
optimum dimensions of circular fi^s of trapezoidal profile 
with variable thermal conductivity and heat transfer 
coefficient . Because of the complex nature of the 
governing differential equation, most of the solutions are 
numerical. 

1.3 SCOPS OB THE PRESENT WOBK 

The literature review reveals that not much work 
has been done on fins with variable thermal properties. 
Whatever little work has been done, it confines itself to 
the study on fin performance of only one thermal property 
at a time, either the variable thermal conductivity or the 
variable heat transfer coefficient. Proper knowledge of the 
effect of each of these propertieg separately on fin 
performance is, no doubt, quite essential for accurate design 
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purposes. Since in actual practice both thermal conductivity 
and heat transfer coefficient come into play simultaneously, 
it is but natural to investigate the combined effect of 
these two variables on fin efficiency. 

Also in most of the papers cited above temperature 
dependence of thermal conductivity has been assumed to 
be linear, A survey of the literature on thermal conductivity 
of metals reveals that a parabolic variation is not only 
more in order but it also accurately reflects the nature of 
the thermal conductivity of metals over a large range of 
temperatures. Keeping this in view, analyses have been 
made to predict both the individual and the combined 
effects of these thermal properties on fin-solutions . A 
parabolic variation for the thermal conductivity and a 
power series function for the variable heat transfer 
coefficient has been assumed. 
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CHAPTER 2 

PRQELSM RORI4I31ATION AW SORUTION 

2.1 STATEMENT OP TH3 PROBLEM ARP THE ASSUPIPTIORS IHVOLVEP 

Accurate knowledge of temperature variation and 
heat flux along the length of a fin is imperative for its 
proper design. The governing heat transfer differential 
equation for a fin has been solved by many authors with 
certain simplifying assumptions. The work presented herein 
entails the solution of the governing equation with some of 
the assLEiptions rd,axed and the comparison of results with 
the conventional solution available from literature. 

The conventional mathematical analyses are normally 
based upon the following assumptions credited to Murray [1] 
and Gardner [2] , 

1. Heat flow is steady 

2. Pin material is homogeneous and isotropic 

3* There are no heat sources or sinks in the fin its^f 

4. Heat flow to or from the fin surface at any 
point is directly proportional to the temperature 
difference between the surface at that point 

and the sihrr.qunding fluid 

5, Thermal cord Uiptiy^'^y of the fin material is 

constant \ , ' ■ ' • 



6 


. Heat transfer coefficient is the same all OTer 
the fin surface 

7, Temperature of the surrounding fluid is uniform 

8, Temperature at the base of the fin is uniform 

9, J’in thickness is so small compared to its height 
that temperature gradients normal to the surface 
may be ne^ected 

10. Heat transferred through the outermost edge of 
the fin is negligible compared to that passing 
through the sides. 

2.2 RBUIXATION OH SOME OH THB ASSTMPTIONS AM) DERIYiTIOM 
OH liK^UATIOl 

In the present analysis an attempt has been made 
to relax assumptions five and six given above. In forced 
convection the heat transfer coefficient is seen to vary 
all along the fin length. Due to complexities in fin- 
geometry, it becomes very difficult to predict the nature 
of this variation. However, attempts made in this direction 
[ 4-7 ] have shown its effect on the overall heat transfer 
from a fin to be significant. But so far no attempt 
has been made to study the simultaneous effect of variable 
theriaal cormluctivity and heat transfer coefficient on 
the p erf cr mane e of a fin. It is. common knowledge that 
thermal conductivity varies not only with the material 
but also with temperature. Mathematically, variable heat 



transfer coefficient and thermal conductivity can be 
expressed as 

h = h(x), 

and 

K = K (e) . 

Next the general differential equation governing 
heat transfer in fins relaxing assumptions five, and 
six , will be derived. 

The differential equation for the fin' profile shown in 
Fig .-(2,1) is formulated from a consideration of the steady- 
state heat balance over a differential element of length dx. 
The difference between the heat entering the differential 
element by conduction at x ard that leaving the differential 
element by conduction at x + dx is given by 

aQ = [K(e) A(x) g]dx 

To satisfy the heat balance equation, an equal 
amount of heat should be dissipated by convection to 
the surrounding medium. Thus, 

dQ = h(x) (T - T,)dS 

CL 

where dS is the surface area of the fin between x and 
X + dx. Also, 

dS = P(x:) dx which ' reduces to 
dS = P dx 


( 2 , 1 ) 


( 2 , 2 ) 



for a constant cross-sectional area fin. Here P is a 


constant unlike P(x) which is a function of x. 

Prom Bqs, (2, 1 ) and (2,2) 

^ [K(e) A(x) dx - h(x) (T-T^) dS = 0 (2*5) 

or, 




— \ — / 


or, 


K 


(e)A(x) ^+K(e)||a^ 


dx‘ 


dx dx 

- b(^) e § = 0 


or, 


d^e 

dx^ 


1 ^ dK(e) 1 ^ dA(x) 

K(6; dx dx a(x') dx dx 


h(x) M A _ n 

K(e) A(x) dx ® ~ ^ 


In order to non-dimensionalize the above equation. 


we put 


and 


e 


9 / 00 . 


X = x/L. 
This gives 


2 ' 
d e 
— 2* 
dx» 


'> dK(9?) , 1 de* dA(x*) 


mi 


dx dx 


AM 


dx dx 




kTTOT?) S 


e s= 0 


(2.4) 
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Since our primary interest lies in studying how 
stron^y the variable thermal properties affect the fin 
solution, the general term of variable cross-sectional 
area in the above equation could be dropped, lq,(2,4-) then 
reduces to 


d^e . 1 de dK(e) 

? ttTaT — ~ — » 

dx’'^ dx 




= 0 


( 2 . 5 ) ■ 


The variation of thermal conductivity of metals with 
temperature is highly irregular and different for different 
materials, For most of the materials, the variation is 
more pronounced at low temperatures, and becomes moderate 
around room tanperatures and above, as in the case of iron 
and its alloys, or even ne^igible as in the case of aluminium. 
Quite cumbersome expressions for thermal conductivity, 
varying in nature from material to material, are available 
for quite a few materials, which describe its behaviour 
over a large range of temperatures. But fortunately the 
behaviour of some of the widely used fin materials like 
copper, aluminium and iron is smooth at and above room 
temperature and can be very well represented by the general 
parabolic curve given by the following expression 

K(e) = (1 +ae+p e^) (2.6) 

where is the thermal conductivity at room temperature and 
a and p are- the temperature coefficients of thermal 
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conductivity. In the light of the above, it becomes 
imperative to use the general expression for the theimal 
conductivity given by Eq, (2,6) as against holding it constant. 
Constant thermal conductivity assumption may not only 
introduce significant errors into the solution but also 
restrict its scope, Sq. (2.6) can also be written as 


K(e) = Kq (1 +«®o ^ *^®o ^ 

or 

K(e’ ) = Kq( 1 + a' e-' + (2.7) 

where a’ and P' are the non-dimensional temperature 
coefficients of thermal conductivity 


and 


a' = “ 


«' 2 
P = |3 . 


Substitution of Eq. (2 .7) inEq.(2,5) yields 


2 

d 9» ^ 

-h 

dx' 


K d+a' 9'+ p' 9'^) dx' ° 


K (1 + a' ©' + p' e'^) 


0. = 0 

K^d + a' 9' + P' 9' )A 


or, 
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d^e» , .des 2 , 2 p' e' /de \2 

— ^ + 2 V — T> + i — f V — r; 

dx' (l+a'^'+P'e’) dx (1+a©+p’0-') dx 

_ ^ 0 
K^A (1+a> e» + P»e'^) 

( 2 . 8 ) 

If a’ and p ' are put equal to zero in the above 
equation, it reduces to the simple case of constant thermal 
conductivity, viz, 

e. = 0 (2.9) 

dx» % ^ 

Comparison of Eqs. (2.8) and (2.9) reveals the 
presence of two additional terms inSq, (2,8) viz. the second 
and the third which take into account the effect of temper- 
ature dependence of thermal conductivity. In order to 
ascertain the significance of these terms, an order of 
magnitude analysis has been resorted to. The analysis 
requires a fair knowledge of the approximate values of a' 
and p' in the range of temperatures involved. Table 2,1 
below gives the approximate values of thermal conductivity 
at different temperatures for copper, iron , and austenitic 
steel [ 15 ] . 
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TABLE 2.1 

VALIBS OB THH1RI4AL GQEDUGTIV ITY AT LIBPERElSrT 
TEI-CPERAfURBS IN 


Temperature 

(®K) 


Material 



Copper 

370 

410 

400 

390. 

370 

Iron 

100 

85 

75 

65 

55 

Austenitic Steel 

9.5 

12,0 

14.5 

17.0 

20.0 


Curve fitting in the range of 300°-500°K yielded the 
following values of K , and A for the above three 
materials (see Table z.2), ' 


TABLE 2.2 

VALUES OB a,p AML BOR DIBBBREMT MATERIALS 


r 

! 

Materials { 

t 

1 * 

K (W/m^K) i 

0 f 

a(l/oK) 

r 

i P(i/°K^) 

_J 

Copper 

390 

-1/8000 

-1/80,0000 

Iron 

75 

-1/750 

0 

Austenitic 

Stainless 

Steel 

14.5 


1 

60,0000 



13 


If the maximum value of a aj^e is taken to be 1/67O 
then a' canes out to be on the order of 0.3 for a temper- 
ature excess of 150®K. A purely parabolic temperature 
dependence of thermal conductivity would require p ' also 
to be equal to 0,3 for the 30 percent increase in the value 
of thermal conductivity. The maximum value of 0-’ being 1 , 
the typical values of coefficients for the second and 
the third terms viz, a‘/(1 + <x' O' + P’ ©■* ) and 
2p'0'/(1 +a» O' +|3'0'^) cornea#^ out to be 0,23 and 0.46 
respectively. Since the coefficient of the first term is 1, 
and the second and the third terms are also quite significant 
as computed above, they can not be ignored and must therefore 
be considered in the analysis. 

YARIAKLB HBAT TRANSFER COEPPICIENT 

literature on the behaviour of fins exposed to steady , 
state surroundings with a variable heat transfer coefficient 
is meagre. The longitudinal fin of rectangular profile 
was studied by Han and lefkowitz [4 ] and Chen and 
Zyskowski[5 ]. Experimental • investigations of Ghai and 
Jacob L6J aM Stachiewicz/suggest a power-series form of 
heat transfer coefficient for rectilinear fins. These 
papers confirm that for practical applications the assumption 
of a constant heat transfer coefficient is unrealistic and 
leads to erroneous results, Han and Lefkowitz recommend a 
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power-series form of h(x) such as 

h(x) = (y + 1) \ (x’ f (2.10) 

where h^^ is the average heat transfer coefficient with 
respect to the length of the fin. (It is to he noted here 
that average value of heat transfer coefficient with respect 
to heat dissipation, h^g^, is different from h^^, in this 
case hgq^ turns out to he greater than hg^,. For details please 
see Appendix C), The reccranended values of v hy Han ard 
Lefkowitz are one and two, which imply a linear and a 
parabolic variation of the heat transfer coefficient with 
the fin length respectively. 

Chen and Zyskowski [5] assumed an exponential variation 
of the heat transfer coefficient given hy 


h(x) 


h. 


1 - a e 


-e X 


a 


1 - a/c (1 - e °) 

Since the experimental results [6-?] confirm the 


( 2.11 ) 


power series form of a heat transfer coefficient given hy 
Fq. (2.10), it is adopted here to study its impact on 
fin performance. Substitution of Bq. (2,10) into Bq, (2,8) 
yields 
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a^6' ^ ^ i—f + 26' »' ^a&' j2 

dx'^ (l + ^e'+P' 0-'^) (I-HZ'&'+P' 0'^) 


l‘ (v+ 1) h p 

( —) 

Ko ^ ■ 


9 "’ X' 


(1 +a’ 9’+^» 9* ) 


=a 0 


( 2 . 12 ) 


Boundary conditions; 9' (0) = 1j d9’{1)/dx' = 0 

where 0< 9* < 1 ; 0 <x’ < 1 

Now the whole problem lies in the solution of ■BqL»(2,12) 
with the specified boundary coniitions.. 

2.3 SOLUTION BOR RECTANGULAR PIN 

Variation of the heat transfer coefficient has been 
assumed as 

h(x') = hg^(2^+ 1) (x')^ 

Bor y = 0 this expression reduces to 

h(x) = hg^ 

and describes the constant heat transfer coefficient case. 

In order to make meaningful comparisons of results with 
the constant heat transfer coefficient case, average value 
of h(x) over the fin length should be found out. Now 

1 1 

^av “ h(x)dxV / dx' = \ 
o o 


Bor rectangular fins 
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A = 2 aM 

P = 2W + 4-7^ = 2W (1 + 2yyW) 

cs>. 2W since 2y^/W « 1 . 


Let, 


2 

L h P 
^ a 


B 


(2.15) 


M' f 

where B = L ^ - y -7 is a non-dimensional number like Biot 
Humber (= h^L/K^) and like it is proportional to the ratio 
of the internal resistance to heat flow by the fin material 
to the film resistance to heat transfer at the surface of 
the fin. The application of fins for heat dissipation 
in any system loses its meaning when Biot Number is 
greater than one because it results in lesser heat dissi- 
pation. But since B = \r(Biot No. x l/y^) and L/y^ is 
usually yery large the upper limit of one does not hold 


good for the parameter B..B cannot be infinitely large 
because other things being equal, the efficiency of fins ■ 
goes down with increase in its magnitude and it, therefore, 
becomes logical to impose an upper limit on its values. 

In the present analysis the upper limit for B was fixed at 
4 based upon the fin dimensions given by Kraus and Snider 
[16]. 


Substitution of B in Sq. (2.12) gives 
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+ a* (d0’/dx' )V(1+a'&'+ »e'^)+2p'e' (d0‘/dx' )V(1+a'^’+P'^’^ ) 

dx’^ 

2 , 0 

~(^+i)b X e'/(i +ate»+ p’e* ) = o (2,14) 

Boundary conditions; 0' (O) = 1; d0' (1 )/dx< = 0 . 

This is a boundary value non-linear differential 
equation of the second order and because of its complex 
nature its analytical scd-Ution is hither to unknown. Bower 
series method of solution also breaks down in this case as 
the resulting series turns out to be diverging in nature. 

Under the circumstances numerical solution remains the only 
way out, A considerable amount of literature is available on 
the numerical solution of differential equations. The 
details of the technique employed to solve the above equation 
and its justification as also the accuracy of the results 
obtained are discussed in Appendix -A. The solution ccxnprises 
the values of 9’ and d0'/<ix' s-t the grid points along the 
fin length, and the fin efficiency for a given set of 
parameters. The differential equation involves 4 parameters 
viz, a', PS ^ ® each set of values for these 

parameters results in a unique solution. The values of the 
parameters with their limits have been listed in Table 2.3 
and the solutions have been obtained for all possiliLe 
combination of the values therein. 



TABLE 2.3 


VALUES OF PARAMETERS 


- 0 . 3 , - 0 . 2 , - 0 . 1 , 0 , 0 , 1 , 0 , 2 , 0 . 5 , 


- 0 . 3 , - 0 . 2 , - 0 . 1 , 0 , 0 . 1 , 0 . 2 , 0 . 3 ., 


V I 0,1,2 


B I 0, 0.5, 1.0, 1.5, 2.0, 2.5, 3, 3.5, 4.0 


For a fuller understanding of the fin and its 
performance, the concept of efficiency comes handy. It- is 
defined as follows j 


Actual heat dissipation from a fin 

Heat dissipation from the fin if its entire 
surface is maintained at the fin root 
temperature 



(2.15 


The heat dissipated from the fin is equal to heat flux at 
the fin-root. Erom Eourier*s Law 



0 (2.16) 


Q 


= - K A 


§1 

dx 


= “ K A 


M = - 

dx 


KAe^ 

L 


de'/dx' 


at X = 


-’^a) S (2.17) 


where hg^ is the average value of the heat transfer coefficient, 
from Eqs.( 2.15), (2.16) and (2.17) 


K A e. 


ri 


- (de'/dx') = - =rT^ (de'/dx') at X = 0 


L h^se 

a o 


K A 


for rectangular fins; 

(Estfer I'lg.2.2(a)) ® 


I h^S 


K ay^w 

1 h 2 IW 

B, 


(2.18) 


y 


(1 +a>e’ +p’ e’^)/B^ 


(2.19) 


from Eqs. (2,18) and (2.19) 

2 2 

^rect “(1+«'6’’ +P'©'' ^ (de’/dx' )/B at x = 0 

= -(1+a> +p' )(de' (0)/dx' )/B^ (2.20) 

(Since 0' (O) = 1 ) 

Also, B,a', P’ and de'/dx' at x = 0 being known, 
the can be readily computed. 
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2.4 SOLUTION gpR SPIITBS (Refer Pig.2.2(b )) 


For a oircular spine of constant cross sectional 


area 


L^h(x)P ^ ^ 2 l\(2^1)x'^ 

KQ(1+a’0' + K^y^j d+a'©-’ + p'©' ) 

(2.21) 

Eqs.(2.12) and (2.21) give 

d^e’/dx'^ + a' (de’/dx' )V(1+a’©'+p'©'^) 

+ 2fue' (d©*/dx’ )^/(1 +a’&’+p’©'^) 

-2B^ (v+ 1) x'^e«/(1 +5^'©' +p’©'^) =0 (2,22) 

Boundary conditions", ©■'(©) = 1j d0''C*)/dx' = 0. 

And, 

^ ^ K^(1+a'e» +p'0«^) y^ Kj^ (l + a’G-' + p'O’^) 

. tot "" oT ^ 

a h„ 1 2 L y_^ 2 ha L 

cL O 

= (1 +«' ©• + p '©’^ )/2B^ 

Substitution of the above expression in Eq, (2,18) 


results in 



ri . = - (1+a' e’ +p' 0-'^ ) (de'/dx' )/2B^ at x = 0 

spni0 ^ 

= -(1 +a' +p’ )(de’jfO)/dx)/ 2B^ ( 2 . 23 ) 

When Bq, (2. 14) is compared with Bq. (2,22 ) ; and Bq,(2,20) 

with Bq. (2,23 ) it becomes obvious that they are similar except 
2 2 

that 2B appears in place of B in case of spines. This leads 
to the conclusion that the solution for both the cases will 
remain the same, except that the scale of B will be modified to 
B/f'Z for spines. The same type of anslysis can be extended to 
deal with the cases where cross-sections of spines are 
regular -polygons. All such cases involve just a change of 
scales, the magnitude of which depends upon the number of 
sides of the polygon. 


'i ^ P n » 





CHAPTER 3 


RESULTS AND DISCUSSION 


In the wake of the failure of analytical methods, 
numerical technique has been eTi 5 )loyed to solve the two- 
point boundary-value differential equation, Eourth Order 
Eunga-Kutta method coupled with Shooting method has been 
preferred to other methods on account of its simplicity 
and accuracy. The division of non-dimensional length 
of the fin has been restricted to 10 grids only because 
further increasein the number of grids adds to accuracy 
oriLy marginally. The values of non-dimensional temperature 
and its gradient ’■'ave been determined at each grid point along 
with the efficiency of the fin using the expression given 
by Eq. (2.20) for a given set of parameters. Solutions have 
been obtained using the simplified assumption of the 
insulated fin end due to the following reasons! 

1 , In no way does the above assumption impede in our 
objective of studying the effect of variable 
thermal properties, 

2, Inaccuracies introduced are insignificant, 

3, It not only simplifies the boundary conditions 
but also eliminates Biot No, (= hg^l/KQ) parameter 
which, if present, would restrict the generality 

of the solution. The Biot No. arises from the real 





boundary condition viz, -Kde/dx = h© at x=Ii 
and it vanishes if the fin end is assumed to 
be insulated. 

The iORTRAN program used to solve the equation satisfies 
the insulated end condition (d©’/dx’ = 0) upto three places 
of decimal and can suitably be modified for further accuracies 
if the situation so warrants. To facilitate a better 
understanding of the results, it becomes necessary to group 
the variables appearing in the Bq.(2.14) and the subsequent 
analyses under the following two heads; 

1. INRSPSNDBIIT VARIABLES; (i) x> , (ii) B, (iii)^', 

(iv) p ’ , and (r) v , 

2. VARIAHI.es ; (i) 6', (ii) de'/dx», 

and (iii) R • 

All the graphical results presented herein have been 
grouped into three categorjes, viz. 

1 , vs. X' with a’, P ‘ ^ ® parameters 

2, d©V^3c’ vs, X' with a Vj and B as parameters' 

3, -I, vs B with a', p', and v as parameters. 

Results falling under category 1 and 2 have been 
drawn on the same graph sheets to enable a clearer presenta- 
tion and a better elucidation. This is so because ©' and 

d©’/dx‘ are interlinked and nature of one can be readily 
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deduced frc*n that of the other. Since the efficiency 
of the fin involves temperature gradient and the temperature- 
modified conductivity, a discussion of temperature and its 
gradient first will lead to a clear insight into the changes 
induced in the s dution. 

The results of the present investigations are presented 

as plots between various parameters from figs, 3,1 to 3.7.' 

* I 0 

figs, 3. 1(a) and 3.2(b) show the variation of non-dimensional 
temperature gradient and temperature with non-dimensional fin 
length for constant heat transfer coefficient condition 
and B = 1 , The curves also demonstrate the effect of the 
temperature dependence of thermal conductivity. The increase 
in the tip-temperature is observed to be on the order of 
17-18 percent for the condition a’ =p’ = 0.3 as against 
the constant thermal conductivity case, a = p =0, This 
is only as expected because positive change in the values of 
parameters a’ and p’ result naturally in an increase in 
the average value of K which further tends to make the 
temperature distribution along the fin more or less uniform. 
The effect of variation of temperature is seen to 

be significant compared to the corresponding effect of 
The reason for this lies in the nature of the rdation c£ 
the variation of thermal conductivity with temperature viz, , 

K(e) = K^(1+ a'e» + p’e‘^)| o<e' < i 



where S'' is the non"-d.itiiensional temperature. Because of 
the magnitude of the 0-* in the above equation, 0-' is 
always less than 0* and thus K increases more rapidly 
with a' than p'. In I'ig.3.1(a) it is interesting to note 
that the tanperature gradient at any given point decreases 
with increase in values of a' and p ' . The decrease at 
the fin-root is on the order of 30 percent for a‘ = P’ = 0,3. 
This result is in consonance with the temperature distribution 
along the fin length. The higher values of thermal conducti- 
vity tend to bring about a unifoim temperature distribution 
along the fin length which in turn results in a decrease of 
the temperature gradient at the fin-root. The small 
difference in the variation of temperature gradient observed 
in the curves f cr identiceO. changes in a' and p' is also 
as expected and can be explained as before due to the 
variation of 0’ fron 0 to 1 . 

The variation of d0'/dx' and x' can be further 
elucidated by considering the simple case of constant h 
and K. The expression for the temperature gradient in 
this case is (See Appendix B) 


dx' 


B tanh B v’-^ere B 


Lv 


K 7 
0 *^ o 


In this equation both B and tanh B decrease 
with any increase in the vilue of K and. since the magnitude 
of K decreases along the 3 cngth of the fin the temperature 
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gradient registers a maximum drop at the fin root. 

Figures 3.2(a) and (b) give the variation of non- 
dimensional temperature gradient and temperature with the 
length of the fin for both positive and negative values of a' 
for constant heat transfer coefficient conditions, B = 1, 
and p= 0, respectively. The temperature gradient is seen 
to increase while the temperature is observed to fall with 
iR03?e€ts^e in 1- ongt h. The effect of negative values of a ' 
is much more pronounced in both the figures. This is as 
expected looking at the nature of the differential equation 
(2J2). 

Figures 3.3(a) and 3.3(b) depict the variation of S'* 
and d9-'/dx* for various values of a', linear variation of 
heat transfer coefficient, and p = 0, The temperature 
gradient is seen to decrease at the fin-root by approxi- 
mately .’16-17. percent. The nature of the ten^erature 
gradient curve in this case is convex upward as against 
the concave upward shape of the curves for the simple 
case of constant heat transfer coefficient and thermal 
conductivity. G-radient changes from almost zero at the 
root to some finite value at the fin-tip. Since linear 
variation of h necessitates a lower heat dissipation near 
the fin-root and a higher heat loss towards the fin— tip, 
it explains the initially flat nature of the temperature 
gradient. 



The marked drop in temperature gradient at the fin 
root can he explained hy the fact that the heat transfer 
coefficient is zero at the fin-root and close to zero at 
the points near the fin root where the temperature is 
the highest. This results in smaller value of equivalent 
heat transfer coefficient (For details please see Appendix C) 

To explain the sharp drop in temperature, heat transfer 
to and from the fin must be considered. Lower heat transfer 
to the fin, because of lower efficiency, results in higher 
temperature drop in this case whereas lower heat dissipation 
results in a correspondin^y lower temperature drop. Both 
these opposite factors namely lower heat transfer into 
the fin and lower heat dissipation from the fin, are 
operating in this case on account of the smaller equivalent 
value of h. But the effect of the former, i.e. the lower 
heat transfer into tho fin, seems to predominate thereby 
resulting in a higher temperature drop. 

Figures 3.4(a) and 3.4(b) demonstrate the effect of v 
on the temperature gfadiont and temperature distribution 
of a fin- fcxr B=1 and constant thermal conductivity conditions 
With the increase in the value of v both the temperature 
gradient and the temperature sxoeBS • a3?e seen to 
decroase along the fin. Changes seem to be more pronounced 
in the range of v varying from 0 to 1 . The, temperature 
distribution curve seems to undergo a small change upto 



about X' =0,4 and thereafter difference for different 
values of y seems to increase, which becomes maximum near 
the fin-tip. Similarly for the values of the temperature 
gradient, the slope near the fin-root seems to be zero for 
the values of v - ^ and 2 and thereafter it gradually 
increases towards the end of the fin. This is as expected 
because for values of 2 ^ = 1 and 2 the heat dissipation at 
the root would be less in that order. 

Figures 3.5(a) and 3.5 (b) exhibit clearly the 
combined effect of variable thermal conductivity and heat 
transfer coefficient on the temperature gradient and 
temperature distribution variation along the length of 
the fin for B = 1. In Pig, 3. 5(b) the effects of increase 
in the values of K and h seem to work at cross purposes 
with the result that the tip temperature seems to be 
the same for the conditions of ' =p' = .3, v = 2 and 

pi _o, V = 0, Otherwise the difference over 
the ranaining length of the fin seems to be quite significant 
leading to the conclusion that higher values of v and a 
and p load to lower dissipation of heat from the fin. 

Figures 3.6(a) through to 3.8(b) show the effect of B 
on the temperature gradient and temperature distribution 
variation along the fin length for constant thermal conductivity, 
constant heat transfer coefficient and variable heat transfer 
coefficient conditions, both linear and parabdLic, B seems to 



have quite a significant influence on both the temperature 
gradient and the temperature distribution. In all the 
three cases the end temperature seems to drop quite signi- 
ficantly. In Pig, 3, 6(b) the temperature excess at the 
fin tip is only 3 percent of that at the fin root when B = 4-. 
Since large B values are equivalent to higher values of 
heat transfer coefficient, they naturally result in larger 
heat dissipation. This is demonstrated, clearly in Pig, 3. 6 (a) 
by the higher values of temperature^at the fin root, lesser 
heat dissipation on account of smaller temperature excess 
towards the tip of the fin explains the change in the 
slope of the temperature gradient curve frcxn higher to 
lower values along the length of the fin. Por the case of 
y = 1 (Pigs, 3, 7(a) , 3.7(b)) the temperature distribution 

curves follow the same pattern as in the previous case i*e, 
y = 0 but the shapes of the curves representing temperature 

gradient are somewhat different being convex upward initially 
and changing gradually to concave upwards. In this case 
the heat transfer coefficient increases from a minimum 
value at the root of the fin to a maximum value at 
the fin-tip while G-' decreases from a maximum at the fin- 
root to a minimum at the fin tip. The maximum of the 
product h G’ which is proportional to the heat dissipation 
from the fin could lie anywhere in between the fin root 
and the tip. Since the heat dissipation is responsible fccr 
the change in the slope of the temperature gradient curves, 
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^ qeems to bold the key to the unusual 

the above explanation seem 

. They seem to have the maxima of 

behaviour of the curve. J-n o' 

their elopee somewhere in the middle between the root and 
the tip. It is also interesting to note that for v = 2, 
the maocimum value for the slope occurs at some Intermediate 
point. Since heat flux at any point is given by 

^2. de 

Q TT ^ =-K (1+aO +§ & ) 


dx 


dx 


Tr,ncr the fin length because of continuous 
As Q decreases along tne iixi & 

u, ^ +-hP fin to the surroundings, so do 
heat dissipation frcm the 

9 and K also. However, if K decreases at a faster rate 

than Q, If Will obviously tend to increase in accordance 

-p s.bo’V^Q x^0l3**bion* 

with the dictates of tho 

i. x 1 A. ffive the variation of efficiency 
Figures 3.9 to 3.^4 

cinps of aS 6’ ^ efficiency 

with B for different values p 

Which IS a strong function of the parameter B. being cen 
percent at B = 0. decreases asymptoticaily to ssro values 

as B tends to infinity* 

^ the effect of variable conductivity 

Figure 3.9 shows the 

..p. ■ percentage change in efficiency increases 

on efficiency, porceni^ & 

n values of B but flatens towards 
very rapidly for small vaiu 

u, -D For the condition of a* - 
the higher values of -d* 

nit starts from 0 at B = 0 and reaches 
p« = 0 and r « 0 ® 
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approximately 8,5 and 9.5 percent values for B = 2 and B = 4 
respectively, while for a ' = 0, j3 ' = 0.3 and = 0, 

the percentage change in efficiency assumes a value of 6.4 ord-y 
at B = 4, thereby establishing the dominance of over p’. 
This is in consonance with the pattern of results obtained 
aM described previously, 

Figure 3.10 brings out the relative importance of 
negative and positive coefficients of thermal conductivity. 

The negative temperature coefficient of thermal conductivity 
is seen to exercise a more pronounced effect on the fin 
efficiency in the negative direction compared to the positive 
temperature coefficient of thermal conductivity. For a’= 0.2 
and B = 4 the percentage change in efficiency is on the order 
of 6.5 while it is approximately 10.0 for a*= -0.2 and B = 4. 
This can be explained as follows. 

For p ' = 0 and ©’ = 1 (non-dimensional temperature 
at the tip) the coefficient of the second term in differential 
Bq, (2,14) reduces to a‘/('^+a'). Its values are 0.167 and 0.25 
for a* = 0.2 and -0,2 re^ actively. Since the magnitude 
of this coefficient viz, a'/(1+a') represents the strength 

of the temperature <^®P'"^®^e|a.tive coelf icieS^^^ 

the main differential equation,/ is bound to exercise a 

predominant influence on the general solution and the efficiency 

of the fin. 
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Figure 3.11 shows the effect of a ' on efficiency for 
p • = 0,2 and 3^= 0, For a’ = 0,3 and B = 4 the percentage 

change in efficiency is on the order of 8,5 while it was 

9.5 -for a' = 0.3, p' = 0 as seen in Fig, (3,9). This shows 
that any increase in the value of p' has the influence of 
reducing the effect a’ on efficiency and vice-versa. 

Figure 3.12 demonstrates the effect of a * on fin 
efficiency with v= 1, and p » = 0, The percentage change 
in T) for the situation of a* = 0,3, B = 4,0 and v = ^ 
is 12,5 while it is only 9.5 for a’ = 0,3, B = 4.0, and 
V = 0 as seem previously in Fig. (3.9). This shows that 

with variable heat transfer coefficient, the temperature 
dependent conductivity affects the performance of the fin 
more stron^y and positively. 

Figure 3.13 presents the effect of variable conductivity 
and heat transfer coefficient on fin performance. Holding 
the thermal conductivity constant a linear change of the 
heat transfer coefficient brings about a drop of approximately 
40 porcont in the efficiency, while a parabolic variation 
of h results in a higher drop of approximately 56 percent 
at B = 4. Again at B = 4, holding the heat transfer 
coefficient constant and varying the thermal conductivity, 
the augmentation in efficiency is seen to be only about 

15.5 peroer*. This clearly shows the. predominant effect of 
variable heat transfer coefficient on fin perfpimance. If 
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both K and h axe changed simultaneously, the figures for 
percentage change in the efficiency are seen to decrease 
bo 30 and 45 for 3 ; = 1 and v - 2 respectively. This shows 
that both the variables h and X play a strong role in 
affecting the fin performance and none of them should 
therefore be ne^ected. 

Lastly Fig, 3 . 14 compares the efficiency of a rectangular 
fin and a circular spine. As explained before, the scale 
of B changes to for a spine, lower efficiency of 

spines for a given value of B is due to the fact that for a 
given cross sectional area and volume, the surface area of 
a cylinder is minimum. This results in lesser heat dissi- 
pation and hence smaller efficiency for the spines. 



CHAPTER 4 


CONCLUSIONS 


Prom the results of the present analysis, the 
following inferences can be drawn: 

(1 ) The temperature dependence of thermal conductivity, 
though being very small for Aluminium, is quite 
significant in the case of Iron, Copper and 
Austenitic Stainless Steel. Iron and copper 
being widely used fin materials, the study of 
variable thermal conductivity with its resultant 
impact on the perfoimance of fins becomes 
imperative, 

* (2) Both increase and decrease of thermal conductivity 
of metals with temperature occur in practice 
depeniing upon the material and the range of 
temperatures involved. Temperature dependence of 
thermal conductivity around and above room 
temperature, for most of the materials, can be 
expressed by the following parabolic expression 

K(0-) = Kq (1 + a0- + 

(3) The rise in thermal conductivity with temperature 
results in 

(a) a more uniform temperature distribution 
along the fin length 
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(b) a drop in the temperature gradient at the 
fin-root 

(c) an increase in the efficiency of the fin* 

The reverse of (a), (b), and (c) happens in case of a drop 
in the value of thermal conductivity with temperature. 

(4) The p)ercentage change brought about in the 
efficiency of a fin by a certain change in 
the value of its thermal conductivity is not 
proportional but less. For example, a 30 percent 
change in K results in an efficiency change of 
only 10 percent. iU. though less it is thought 
significant and therefore justifies the modified 
analysis. But all this is valid only if the 
parameter B is large (i.e, B^l). For small 
values of B, K can be safely assumed to be constant 
without the risk of introducing any significant 
error into the results. 

(5) Out of a’ and p’ j the former is seen to exercise 
a stronger influence on the fin efficiency. 

(6) The Impact of variable h on fin performance is 
much more pronounced than that of K, and should be 
considered for all possible values of B. The 
assumption of variable h and constant K leads to 
significant errors and it is, therefore, recommended 
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accurate 

that both be treated as variables lor 
results, 

, _T.i o h and K 

(7) The overall combined effect of variao 

-hv as much 

is seen to reduce the fin efficiency oy 
as 65 percent in some cases and so a 
safety factor needs to be employed for pt'oP 

design purposes. 
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CHAPTER 5 

SUGGESTIONS iDR FUTURE NOBK 


The limitation of the present 
applicability only to fins with constant cross-sectional 

area. It also excludes from its purview annular fins 

S3-. y 

on account of their variable cross-sectioh 

present work can be extended on lines as suggested by 

Gardner [2 ] in order to bring under its r^rige all possible 

fin shapes. 

The author failed to pursue the problem along the 
suggested lines of Gardner partly because 

.1. jy o-u rjt'dblBin# How'svgiTj 

and partly du6 to th8 complexxty of tho V 

, area, it is proper 

for the benefit of future workers in thi® 

^ ^ 4 -he method of attack 

to list the difficulties involved and the 


of the problem, 

^ . . ._£, heat transfer in 

The differential equation governing 

fins with variable conductivity, heat transfer 

and cross-sectional area is ^ven as fol^®^®* 


d^e-* 

dx«^ 


+ (1/A 


M .'I it’ 

dx ' dx ' 


^ V' 





^ . M . ’ ±1 ^ (dO’/dx' f 

(l+a'^’+p'e* ) 




(1+a« 8'+P* 8'^) 


0 


(5.1) 



58 


Evaluation of dA/dx involves simple calculus but 

Q Q! 

^ poses a problem, ^ can, however, be calculated 
from the following expression deduced by Gardner [2] 


• h dS 

U = Xv/" rrr- T— 

lA dx 


where u for different cases is given as follows. 


SIRAIGHT PINS: 


Por variation of thickness like 
1 -2n 

Jq 


The value of u is given as 


2(1 -a) 


2. SPINES ~ GIRGIILAR OR POLYGONAL SECTIONS : 

The corresponding quantities in this case are given 


7 ~ 7q (xA) 


1-2n 
i -n 




3. ANNUI.AP PINS: 


y = yo (xA)'l“n 


u = (1-vi) (xA)"^”^ ^ 
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For n = 0, 1/2, 1, 3/?^, and 2 substitution of expressions 


dA dS 

for ^ and ^ for different cases in Bq.(5.1) yields 


+ JT ^ + a* (d9'/dx* )^/(1+a'&'+p«0-«^) 


+2 p'e' (de'/dx* )V(l + a’^‘+ |3‘e‘^) 


x^ (1+a‘ e'+ 13 


(5.2) 


Boundary conditions: (O) = 1; =0, < 

0 <8-‘ < 1 5 0 < x» < 1 , 

where G and D are constantsi n = 0, 1/2, 1, 3/2 and 2, 
h(x) = hg^ f (x) and = B^VVo* 

The presence of second and fifth terms in the above 
equation introduce a signularity at x = 0; and the whole 
problem boils down to one of solving a two -point boundary- 
value non-linear second order differential equation with a 
singularity. The singularity poses a big problem and many 
a numerical techniques breakdown in this case (See Appendix 
A). However, the problem involved is certainly not intract- 
able and its analysis- will go a long way in providing a 
highly generalized solution for fins of all shapes. 

In the absence of experimental data on heat transfer 
coefficient for fins of different geometries, one has 
largely to depend upon theoretical models. A general solution 



of the Eq.(5.2 ) with h(x) based on experimental results 
should be attempted and a Terification of the results 
thus obtained confirmed by subsequent experimental investi- 
gation into the performance of fins. Although a large amount 
of work is involved in the process, it is certainly worthwhile 


to pursue the same. 
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i\PPBiroiX“A 


The differential equation encountered in the present 
analysis is a two-point boundary'-ralue non-linear second 
order equation. There are many numerical techniques to 
handle such problems e.g., Pinite Difference Method, least 
Square Method, Runga-Kutta Method etc. Since Runga-Kutta 
Method handles only the initial value problems, it has been 
coupled with Shooting Method to take care of the boundary- 
value problem. The methoddLogy involves assuming a value 
for Q-'t the fin-root, solving the equation and 

checking whether the boundary condition at the fin-tip 
(dS-*/dx' = 0) is satisfied. If not, then assuming a new value 
for the temperature gradient and thus keep on repeating the 
process till the boundary condition at the fin-tip is 
satisfied. 

The method has been given preference over other ones 

on account of its simplicity and accuracy. The Fourth 

with 

Order Runga-Kutta Method coupled/ the satisfaction of 
the boundary-condition at the fin-tip upto three places of 
decimal gives an excellent accuracy. Table A-1 compares 
exact and numerically obtained values of non-dimensional 
temperature gradient, dO'/dx', at the fin-root for constant 
h and K values. 
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TABLE A-1 




1 

B ! 

1 

t 

f 

f 

0 

1 I 

f t 

1 ^ J 

1 f 

! ! 

f r 

- - t 

2 

f 

f 

I 3 

f 

! 

f 

A .. 



i 4 ’ 

f 

1 

! 

r 

Exact 






d©' 

IF" 

'0 

.76159 

1 .92806 

2.98515 

3.99732 

(A) 






Numeri- 

cally 

obta- 

ined 

0 

.76169 

1 .92852 

2.98549 

3.99746 

d©‘ 

dx^ 






(B) 






Differ- 

ence 

(B-A) 

0.0 

.00010 

0.00046 

0.00034 

0.00014 

Since 

changes of the 

order of 

0.01 to .1 

are involved, 


tlie accuracies obtained are more than sufficient. 




APPBHDIX-B 


SOLUTION POE PINS 


The differential equation governing heat transfer in 
fins with constant heat transfer coefficient and thennal 
conductivity is given as 


where, . 


2 

d O' 

2 

L h P 
a 

2 

dx' 

A K 

0 

d^O« 

dx'^ 

“ O' = 0 

2 

2 

I h P 

CU 

m 

A' K 


Boundary conditions 6'(0) = 1, ^ = 0,. 

Solution of the above equation with its boundary 
conditions yields 

e’ = Cos h m (1 “X' )/Gos h m 


(B.1 ) 


^ “ Cos h m 

Differentiation of Bq. (B. 1 ) gives 

d©'/dx' = -ID Sinh h m (1-x' )/Cos h m 


-m tan h m 
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Now, 


T] = 


K^A ©» (0)/dx' ) A 

P 1? h PITh 

0 a a 


d©' (0) 
dx’ 


1 



d©« (0) 
dx’ 


m tan h m 
2 


m 


= - tanh m/m 


For rectangular fins ; m = B, 
For spines (circular): m = B'/2' , 
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APPENDIX-n 

EQUIVALENT HEAT TRANSEER COEEEICIENT 


The form of the heat transfer ooeffiolent suggested by 
Han and Lefkowitz [4 ] is 

h = (v+ 1)h„ 

d 


It'S average value with respect to the fin length Is 

V = / (’'+ 1) \ = h . 

O O' 

But this value of average does not take Into account the 
tanperature variation along the fin-length. 

It is logical to define equivalent heat transfer 
coefficient as that ooefflciont which, if assumed constant, 
would result in the seme heat dissipation as the variable one. 

Since, Q = h Q(P dx' ). 
o 

We get, 


or, 


1 

/ hg^ & P dx« 


1 

f 

o 


h et p dx' 



1 

/ h 0 dx’ 

0 

1 

j* 0* dx* 
o 


x> 0* dx* 

0 

1 

/ 0’ dx’ 

o 


Por 1^= 1 . 
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1 

2h f X' &‘dx’ 

cL 



sin h m 
m 


1 

2hg^ ^ 6 !os h m(1-x' )x'dx' 
} Gos h m (1-x' )dx» 


2h (Cos h m-l ) 

Q. . 

mSin h m 


The value is 1/2 at m = 0, and 

decreases continuously till it reaches zero at ra - oo , 
h„„. is therefore, less than h^ except when m = 0. 
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